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Abstract 

The introduction of a lattice converts a singular boundary-layer problem in the con- 
tinuum into a regular perturbation problem. However, the continuum limit of the 
discrete problem is extremely nontrivial and is not completely understood. This paper 
examines two singular boundary-layer problems taken from mathematical physics, the 
instanton problem and the Blasius equation, and in each case examines two strate- 
gies, Pade resummation and variational perturbation theory, to recover the solution 
to the continuum problem from the solution to the associated discrete problem. Both 
resummation procedures produce good and interesting results for the two cases, but 
the results still deviate from the exact solutions. To understand the discrepancy a 
comprehensive large-order behavior analysis of the strong-coupling lattice expansions 
for each of the two problems is done. 

I Introduction 

In this paper we report some major advances in understanding (albeit not a complete solu- 
tion to) a difficult general class of problems in mathematical physics. We consider here the 
conversion of a continuum problem into a discrete problem by the insertion of a lattice spac- 
ing parameter a, the solution of the continuum problem on the lattice, and the subsequent 
extremely subtle continuum limit a — > 0. 

Almost every continuum physics problem is singular as a function of the parameters 
in the problem. As a result, only rarely does the perturbation series take the form of a 
Taylor series having a nonzero radius of convergence. As an elementary example, consider 
the algebraic polynomial equation 

ex^ +x-1^0. (1) 
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This problem is singular in the limit e — s- 0. In this limit, the degree of the polynomial 
changes from three to one and thus two of the roots abruptly disappear. As a consequence, 
a perturbative solution to this problem [expressing the roots x{e) as series in powers of e] 
yields expressions that are more complicated than Taylor series. 

A more elaborate example of a singular problem is the time-independent Schrodinger 
equation 

In the classical limit h —^ this differential equation abruptly becomes an algebraic equation, 
and thus the general solution no longer contains any arbitrary constants or functions and, 
as a result, it can no longer satisfy the initial conditions. We know that for small h the 
solution is not Taylor-like but rather is a singular exponential in WKB form: 

*(x) - e^(")/'* (h^O). (3) 

In the study of quantum field theory, it is well known that infinities appear in the 
perturbative expansion in powers of the coupling constant. There are two kinds of infinities. 
The first kind, which is due to the point-like nature of the interaction, requires the use of 
renormalization. The second kind, which is due to singularities in the complex-coupling- 
constant plane, forces the perturbation series to have a zero radius of convergence. 

A quantum field theory can be regulated by introducing a lattice spacing. The resulting 
discrete theory is completely finite and can be studied numerically by using various kinds 
of numerical methods such as Monte Carlo integration. However, the underlying singular 
nature of the continuum quantum field theory resurfaces in the continuum limit a ^ 0. 
The introduction of a lattice spacing and the singular nature of the continuum limit was 
investigated in a series of papers by Bender et al. . 

A quantum field theory is just one instance in which discretization regulates and elimi- 
nates the singular nature of the problem. It is also known that introducing a lattice spac- 
ing converts a boundary-layer problem, which is a singular perturbation problem, into a 
regular perturbation problem [llOH12[. A boundary-layer problem is a differential-equation- 
boundary-value problem in which the highest derivative of the differential equation is mul- 
tiplied by a small parameter e. Consider as an example 

ey"{x) + a{x)y'{x) + b{x)y{x) = c{x), (4) 

where the boundary conditions on the function y{x) typically have a form such as 

y{0) = A, y{l) = B. (5) 

This boundary-value problem is singular because in the limit e ^ one of the solutions 
abruptly disappears and the limiting solution is not able to satisfy the two boundary condi- 
tions in (^. The usual way to solve the boundary- value problem (^) - is to decompose 
the interval < a; < 1 into two regions, an outer region, in which the solution varies slowly 
as a function of x, and an inner region or boundary-layer region, in which the solution varies 
rapidly as a function of x. The boundary-layer region is a narrow region whose thickness is 
typically of order e or some power of e . 

An important example of a boundary-layer problem is the instanton equation 

ey-"{x)+f{x)-f{x)=0, (6) 

with the associated boundary conditions 

/(O) = 0, /(oo) = 1. (7) 
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The exact solution to this instanton problem is 

/(a;) =tanh-^. (8) 

Note that the solution f{x) varies rapidly at the origin a; = over a region of thickness e; 
this is the boundary-layer region. The solution varies slowly (it is approximately 1) outside 
of this region. The outer region consists of those x not near the origin. 

A novel way to solve the instanton problem is to discretize it by introducing a lattice. 
On the lattice, the differential equation becomes a difference equation that can easily be 
solved perturbatively. In the continuum limit, as the lattice spacing vanishes, we then obtain 
a strong-coupling expansion that must be evaluated by means of a Pade or a variational 
perturbation theory method. To illustrate the approach our objective will be to calculate 
the slope of the instanton at a; = 0, which from (^) has the value 

/■ (0) ^ (9) 

We introduce a lattice with lattice spacing a so that the real axis is discretized in steps 
of width a. The spatial coordinate reads Xn = na, where the function f{x) assumes the 
value /„ = f{xn)- On the lattice the second spatial derivative in (^) becomes 

/ [x) 2 ■ 

Thus, from the instanton equation (|^) we obtain the difference equation 

- 2/. + fn-l) + /„ - fn = 0, (11) 

where the boundary values follow from (^: 

/o = 0, /oo = 1. (12) 
The natural expansion parameter now is /a? , to which we assign the name 5: 

The singular perturbation problem in the continuum [whose solution fix) in (^ does not 
possess a Taylor expansion in powers of e] , has become a regular perturbation problem. That 
is, we can now expand the solution /„ to the difference equation ( pT| ) as a Taylor series in 
powers of 5: 

fn = an,o + an.iS + a„^2'5^ + . . . . (14) 
We impose the boundary values ([l^ ) by requiring that 

ao.o = and a„^o = 1 {n > 1). (15) 
Inserting the ansatz (^4[) into the difference equation (pT|), we get the recursion relation ||ic[] 

^ ^ j — l ^ J — 1 j — k 

'^■nj = 2^n+l,j-l + Q^n,i-1 + '^'^n-l.j-l — 0'n,k0.n,j-k — 0,n,kO'n,ian,j -k-l ■ (16) 

fc=l k=l 1=1 
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For the first derivative at the origin a; = this leads to the series 



f'(Q) = lim = lim — ~ lim —} ai 

a-+0 a a^O a a^o a ^ 



We have calculated the coefficients anj with the help of Maple V R7 up to order j = 200. 
The first 20 numbers are given in Table A complete list of these coefficients can be found 



on the webpage of the author FW |14|. Note that the expansion parameter (5 in (17) is not 
small but rather tends to infinity in the limit as the lattice spacing a approaches zero. Using 
the parameter 5 defined in (^3|) we rewrite the series ([l7| ) as 

/■<")^L-^('-|4-w-i as, 

Taking into account the exact result (^), we obtain the identity 

5 11,54 



—= = lim \f5 I 1 — 



V2 <5^oo V 2 8 128 



(19) 



The purpose of this paper is to examine equations like ([191). This equation shows that 
the singular nature of the instanton problem has resurfaced in the continuum limit i5 ^ 00 
of the lattice expansion. The expression on the right side of (^9|) should have the value 
l/-\/2 = 0.7071067812 . . ., but it is not at all obvious why this is so, and the objective of 
this paper is to analyze this difficult and subtle limit. 

This paper is organized as follows. In Sec. |l| we use Fade techniques to perform the 
limit in (p^. We will see that while the results are not bad (the accuracy is about 1%), 
better methods are needed. We perform the Fade analysis to much higher order than 
has ever been done before and we discover a new qualitative behavior that has not yet 



been observed. In Sec. Ill we try the use of the variational perturbation theory techniques 
introduced by Kleinert to perform the sum in (p^). These techniques increase the accuracy 
by a factor of about 10, but they still do not give the exact result. While variational 
perturbation theory works very well in summing the strong-coupling series for the ground- 
state energy of the anharmonic oscillator ]l5[ | , and for the critical exponents of second-order 
phase transitions |]l^, we show that the series in ( p!9| ) is at the very edge of validity for 
Kleinert's methods. We then examine the large-order behavior of the terms of the sum 
( pj| ) in Sec. IV. We show definitively that the Taylor expansion has a nonzero radius 



of convergence and thus, on the lattice, the instanton problem is a regular perturbation 
problem. 

In Sec. we turn to a more difficult singular perturbation problem; namely, the Blasius 
equation of fluid dynamics. We use the same approach as for the instanton equation. In 
Sees. |Vl4 and VIII we study the summation of the lattice perturbation expansion 



using Fade and variational methods and we examine the large-order behavior of the lattice 
perturbation series. We find that Fade methods give good but not excellent results and that 
variational perturbation theory is better than Fade. Again, the series we need to evaluate 
in the continuum limit lies at the very edge of validity for Kleinert's methods. We also find 
that, unlike the lattice perturbation expansion coefficients for the instanton problem, the 
sign pattern of the Blasius weak-coupling series does not alternate. Rather, it is governed 
by a cosine function with a frequency different from tt. 
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II Pade Resummation for the Instanton Equation 

In this section we examine what happens if we attempt to evaluate the right side of ( [l9| ) 
by using Pade techniques. Pade resummation has already been applied to the instanton 
problem up to 50th order ||l0| . However, we have been able to perform the procedures to 
much higher orders. We have discovered that remarkable and unsuspected new phenomena 
occur just a few orders beyond what has been computed. 

The procedure is as follows. Consider the formal Frobenius series 

oo 

S'(5) = ,5^^^a„(5", (20) 

n=0 

where M is a non-negative number. Raising this series to the power 1/Af, inverting the 
right hand side and re-expanding, we obtain 

5i/M(5)^_^, (21) 

with new expansion coefficients 6„ . Assuming we know the first -I- 1 terms of the original 
power series in (|2^) , we raise equation ( |2l] ) to the power N. We then truncate the summation 
at 71 = TV, finally getting 

SN/M^g^ = ^ , (22) 

n=0 

where we have re-expanded and obtained new expansion coefficients c„. In the limit (5 — > cxd, 
only the iVth term in the denominator survives and we obtain the approximant 

iS^f/^' ^ hm S^^'\6) = lim = (23) 



n=0 



The approximant — (^'^jv'^ j zeroth-order survivor of the limiting process. 

Also, taking into account the first-order correction we observe that, as in the case of varia- 
tional perturbation theory (see S ec. |lll|) , there is an approach to scaling. In the limit (5 — > oo 
the Frobenius series S{S) in Eq. (|20|) converges to a constant C. Additionally, the approach 
to scaling, following from the Pade resummation ( p3| ) , reveals how fast it converges: 

S{5) ^ C + C'6-^ {6^oo). (24) 

We now apply this procedure to the boundary-layer problem (ITTI). [Recall that the weak- 
coupling coefficients for the first 20 coefficients aij obtained from (|l6|) are shown in Table 
|l| and that more can be found in jlj].] Rcsumming the series ( |l^ ) for n — 1, 

N 

/i=^ai,,J^ (25) 

3=0 

according to the Pade procedure ( p3| ) with i\/ = 1/2 as follows from ( [l9| ) and evaluating the 
approximants Sn — [cj^ ) , we get the numbers listed in Table g. 
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Compared with the numerical solution l/\/2 w 0.7171067812, this strong-coupling ex- 
pansion seems to converge quite well. However, when we go to higher orders, we find that 
the numbers drop below the exact solution and assume a minimum at A'^ = 24, where the 
approximant has the value S'24 ~ 0.70198319. The approximants then rise again, cross the 
exact solution at iV = 41 and become complex at = 52. The appearance of complex 
numbers is a consequence of taking the A^th root in equation ( p3| ) when the coefficients c^'' 
become negative. This phenomenon has not been observed before in the course of using 
this Padc procedure. The imaginary part then becomes smaller and smaller as N rises. 
Abruptly, at = 68, the approximants become real again. As one can see from the spikes 
in Fig. |l| this pattern is repeated for higher A^. Note that the figure only shows the real part 
of the Fade approximant Sn- 

Apparently, the sequence of approximants Sn does not converge. The singular nature 
of the instanton equation has the effect of making the Fade approximants behave like the 
partial sums of a divergent (asymptotic) series; at first the partial sums appear to converge 
to a limit, and then they veer off. In the case of the Fade's shown in Fig. |l| the approximants 
approach to within 1% of the correct limit before veering off. It appears that another more 
powerful resummation technique is needed to treat the expression in (|l^). In the next 
section we apply a technique due to Kleinert. 



Ill Variational Perturbation Theory for the Instanton 
Equation 

Kleinert has developed a technique in the context of the ground-state energy of the an- 
harmonic oscillator and of critical exponents of second-order phase transitions ||l^ for 
summing divergent perturbation series. This technique, known as Kleinert's square-root 
trick, is described below. 

Consider a weak-coupling series 

N 

which is truncated at order A^. Rewrite this weak-coupling expansion by introducing an 
auxiliary scaling parameter k [|15|,|1^: 

n=0 ^ ' 

which is set to k = 1 later. The square-root trick now reads 



K Vif2 _^ ^2 _ ^2 ^ /^Vl + (28) 
where isT is a "dummy" scaling parameter and 



In the case of the anharmonic oscillator, K is the frequency of a trial harmonic oscillator 
Substituting ( p8| ) into the truncated weak-coupling series (|27|), we obtain 

N 

fN{S, K) = Y^ /„iff-"«(l + fr)(p-"9)/25" . (30) 

n=0 
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The factor {l + gr)°' with a = {p — nq)/2 can be expanded by means of generahzed binomials 
according to 



N-n / s N-n / \ / -, \ fe 



A;=0 ^ ' k=0 ^ ' ^ ' 

where we have used ( p9| ) and finally have set k = 1. The binomial is defined as 

r(a + i) 



^kj r(fc + i)r(a + fc + i) ■ 

We deduce that the function /at (5, K) can now be written as 



N 



E 



n=0 Lk=a 

To first order this expression reduces to 

fiiS,K) = (l - f ) foKP + poKP-^ + fiSKP-^ 
Applying the principle of least sensitivity ||l7|] leaves us with 



dfi{S,K) 
dK 



foK-^ + {p-q)fiSK-'' = 0. 



Next, making the strong-coupling ansatz 

we obtain the following equation from (^5|): 
P{P - 2) 



+ 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



The second term is a subleading contribution in the limit as the coupling d goes to infinity 
which we can neglect. Solving for kg^^ we then get 



(1) _ / 2/i p~q 



V /o P(p-2) 



1/? 



(38) 



Assuming that the ansatz (36) for the variational parameter K{S) also holds for higher 
orders we obtain from the fmiction fpf{S,K) in ( |3^ ) 

MS)^S^ [bi''\klf^^) + briki''\kr)6-y^ + ..] , (39) 
where the leading strong-coupling coefficient bQ^\kjf^'^) is given by 



The inner sum can be further simplified, using 



a — 1 



(40) 



(41) 
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Thus the strong-couphng coefficient (^) reduces to 



N 



c'(fcr)=E(-i) 



N-n ( 2 



n=0 



N-n 



fn{kQ 



p-nq 



(42) 



So, looking at equation ( p9|) we see that the fraction p/g tells us the leading power behavior 
in 5 and 2/q indicates the approach to scaling: 



^ f^5^ ~ dP/t [bo + hS-^/" + ...) (,5 ^ oo). 



(43) 



For the instanton equation we can determine the numbers p and q by re-obtaining the 
differential equation (|6|) from the difference equation (pi]). The positive real axis is dis- 
cretized in steps of width a, so that we let x„ = na. The power series expansion for the 
discrete function /„ = f{xn) has the form 

fn±i = f{xn) ± nx^)a + ^.f"ix„)a^ ± i/"'(x„)a3 + ^f""{x^)a^ ± ... . (44) 

Thus, the numerator of the second derivative ( |l0|) becomes 

1 



/„+i - 2/„ + /„_i = + ^/;"'a4 + ... , 



(45) 



so the zeroth-, first-, and third-order contributions cancel. Translating the lattice result for 
/„ back to the continuous function ,f{xn) = /«, the difference equation ( pi] ) reads 



+ f{x)-f{x)=Q. 



Writing out the power series 

fix) = fo{x) + a^hix) + a^f2ix) + 
and comparing even powers of a, we get from equation ( ^6| ) for a'^ 

e'/o(^) + /o(^)-/o(a^-)=0, 
which is just the original instanton equation (^). For we have 

1 



^Vr(^) + /iW(l-3/o'(x)) 



12 



e'Cix). 



The boundary values read 



/o(0) = 0, /o((X3) = 1, 



(46) 



(47) 



(48) 



(49) 



(50) 



and 



/i(0) = /i(oo) = 0, 

respectively. The solution to equation (El) with the boundary values (^oh is of course 



(51) 



/o(a;) = tanh ■ 



:V2- 



(52) 
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So, finally from 



Comparing equation 
2 



I) we get for the derivative at the origin x — 0: 

8|), we resum the weak-coupling series in (|l8|) as 



/'(o) ^ m) + -m + ... = ^ + -f[{o) 



I with 



1 



52 



-1/2 



1 



(53) 



(54) 



Also, comparing with (| 
are given by 



we conclude that the leading power and the approach to scaling 



P 

q 



1 

2' 



2 

- = 1, 

q 



(55) 



respectively. So we identify p = — 1 and q = 2. 

We now evaluate the leading strong-coupling coefficient bo from ( ^ ) according to ( ^ ) 
with p = — 1 and q — 2. To that end we substitute our 200 weak-coupling coefficients 
from into the formula using a computer algebra program. We are now confronted 
with the following problem: The principle of least sensitivity cannot be unambiguously 
applied. Optimizing with respect to extrema, inflection points, or higher derivatives does 
yield converging results for the strong-coupling limit. However, all these strong-coupling 
series converge to the wrong values. 

There is one particularly unpleasant case: The second derivative with respect to fco for 
the largest kg where this derivative exists (see Fig. ^) gives a convergent strong-coupling 
series. The numbers come extremely close to 1/V2 as one can see from the 20 numbers in 
Table ^. The 200th leading strong-coupling coefficient is 6q^°'''' = 0.707417.... However, a 
Richardson extrapolation based on the first 200 orders then unfortunately shows that 
variational perturbation theory produces a value slightly smaller than 1/V2. The first six 
orders of Richardson extrapolations are presented in Table ^. Hence, the strong-coupling 
series 5q^-* does converge, but it converges to the wrong number, only one part per 1000 
away from the true value: 



/i 



(VPT) 



lini 



200 

E 

n=0 



0.7063998320858845 ± 0.0000000000000001 



(56) 



compared with /'(O) = 1/V2 = 0.7071067812... . The deviation is just 0.099%, but 1/V2 
can unfortunately be ruled out. 

Given that p = —1 and q — 2, the failure of variational perturbation theory is not 
surprising. According to Ref. the fraction 2/q must lie within the open interval (1/2, 1). 
Otherwise, one cannot prove that variational perturbation theory converges. Thus, this 
problem lies exactly on the upper boundary of the region in which the summation method is 
known to work. 

We can understand the upper edge of the range of the parameter 2 /q that describes the 
approach to scaling 2/q by looking at the standard deviation from the actual limiting value. 
It turns out ||l^ that the deviation in the limit as the perturbative order TV goes to infinity 
assumes the shape 



-bn 



bo 



exp 



(57) 



where C is a constant. So, to obtain exponential convergence for the sequence formed by 
the bo^\ we need 1 — 2/g>0. In other words, the approach to scaling 2/q is bounded and 
it must be smaller than one. The lower edge is more subtle and is discussed in Ref. [|6|. 
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In conclusion, we have applied variational perturbation theory to a case that lies at 
the very edge of its applicability. We see that variational perturbation theory gives better 
results by about a factor of 10 than the Pade approximations examined in Sec. ||. However, 
we have not yet found a systematic method for resumming ( [l9| ) that enables us to perform 
the continuum limit of the discrete lattice theory. Therefore, we now lay the foundation for 
further investigations by analyzing the large-order behavior of the instanton series. 



IV Large-Order Behavior for the Instanton Equation 

It can be seen from the numerical results in ]l^ that the instanton weak-coupling series is 
of Borel type. That is, it exhibits an alternating sign pattern. From the ratio test we can 
see that the coefficients a„j- do not grow factorially fast. The large-order behavior of a„j- 
has the general form 

a„,, ^{^ir+^+'Kij^-B^ (j-cx)). (58) 
The constant An can be obtained by evaluating the limit 

g — — 

An = hm (59) 

^(j + i)^ 

and the reciprocal of the radius of convergence is 

i^n - - lim f ^ V" . (60) 

Also, the overall factor i?„ is determined from 

Bn = hm (61) 

Using the 200 weak-coupling coefficients, we find that the exponent An and the reciprocal 
radius of convergence Kn are independent of n. The value of K2 — 2.46682906 coincides 
with Ki = 2.46682906 for aU significant digits. The same is true for Ai = -1.500000 and 
A2 = —1.500000. Thus, it appears that we may omit the subscripts n for Kn and An- 
In contrast, the data suggests that Bn strongly depends on n. En is the numerical value 
associated with the largest uncertainty. In fact, Eq. ( |6l| ) suggests that small deviations in 
K and A lead to dramatic changes in the value of We calculated A, K, Bi, and B2 
up to 200th order with the help of Maple V R7. We then extrapolated these 200 orders to 
infinity using Richardson extrapolation [p3| . We obtained 

A = -1.500000 ±0.000001, 

K = 2.46682906 ±0.0000001, 

Bi = 0.0171 ±0.0001, 

B2 = 0.1190 ±0.0001. (62) 

Detailed numerical results for the first six Richardson extrapolations for the exponent 
A, the inverse radius of convergence K, and the overall factors Bi and B2 can be found 
in Tables ||, |^, 0, and ^. The calculation of Bi is extremely delicate; changing the inverse 
radius of convergence in the sixth decimal place influences the third significant figure of Bi . 
The same is true of B2- 
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Unfortunately, there is no way to derive these values by applying asymptotic analysis 
to the recursion relation ([l^). The problem is that the double summation in this equation 
includes small j, so we cannot let j go to infinity and use the large-order behavior (|5^). 
Substituting the ansatz (|58| ) into equation (|l^) and taking the limit leads to contradictory 
results. For n = I we get 

k=l 

-^^?^EE^-'^^(^'-^-o^- (63) 

k=l 1=1 

Pulling out some factors and letting x = k/j, we obtain for the first summation 

if and only if A > — 1. For A < —1 which is strongly favored by the data we obtain 

dec [a;(l - x)]^ = 2C . (65) 



The double summation reduces to 

-^'—k^i^f^ k i\ _ f\ r\ , , ,^^_r3(A + i) 



k=l 1 = 1 

where y = l/j and A> —1. For A < —1 the result is 



dx / dy[xy{l~x-y)]^ = 3e{~A). (67) 



Substituting the results in (|65|) and ( |67| ) into (|6^) leads to a contradiction: The inverse 
radius of convergence then turns out to be 

1 + 7^ 

""^i + sciDB. + ieiDBr ^''^ 

which would imply that, given Bi = 0.0171 and B2 = 0.1190, the value of K would be 

K = 3.940. (69) 

This result can be ruled out because of the numerical result (|6^). Also, (^) does not 
contain the exponent A because all the factors in ( |6^ ) cancel. So A cannot be determined 
analytically using this asymptotic analysis. 



V Boundary-Layers on the Lattice — Blasius Equation 

The Blasius equation fl^ arises in the study of fluid dynamics. It is a special limiting case 
of the Navier-Stokes equation and determines the flow of an incompressible fluid across a 
semi-infinite flat plate. The equation reads 



2ey"'ix) + yix)y"ix) = 0. 



(70) 
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Assuming that the tangential velocity y'{x) at the outer limit of the boundary layer is 
constant, the boundary conditions read [ |l9| 

2/(0) = y'{0) = 0, y'{^) = 1. (71) 

Our objective here is to calculate the second derivative y"(0), which represents the stress 
on the plate. We discretize the Blasius equation ([70|) by introducing a lattice spacing a: 



2<5(/„+i - 3/„ + 3/„_i - fn-2) + /„(/„+! - 2/„ + U-i) = 0, (72) 
where we define /„ = y{na)/a and S = e/a^. The boundary conditions ( |7l| ) now read 

/o = /-I = 0, /„ {n^ oo). (73) 

Expanding /„ as a series in powers of S as in Eq. ( [1^ ) , we obtain the recursion relation | [Tc[ | 

2 

dn+ij — 2a„j- + a„_ij- — —— (a,i+ij_i — 3a„j_i + 3a„_ij_i — ari-2,j-i) 



— ^ On, A: (fln+l j-fc " 2a„j-_fc + a„_i j_fc) , (74) 



n 
fc=i 

The boundary values are 

a„,o = n (n > 0), 

a-1,0 = 0, 

a-„-ij = a„j (n > 0). (75) 

Eq. (^ can be solved order by order by using a computer algebra program. Table || shows 
the first 20 weak-coupling coefficients aij. All coefficients up to the 300th order can be 
found at Pol. 



VI Pade Resummation for the Blasius Equation 

We now resum the weak-coupling coef ficie nts using the Pade method ( p3|) with M = —1/2. 



This value of M will be derived in Sec. Vll in Eq. (|82|). The exact solution [|10| to the Blasius 
equation (|7^), obtained numerically up to five digits, is y"{0) = 0.33206. Unfortunately, 
the sequence formed by the approximants Sn appears to converge, but not to the correct 
value. According to Table |l^ the sequence becomes very flat and Richardson extrapolation 
shows that the Sn approach the wrong limiting value (see Table |ll|). A third-order 
Richardson gives Soo = 0.3430, based on the first 70 weak-coupling coefRcients. This value 
is significantly higher than the correct value y"(0) = 0.33206, the deviation is 3.3%. 

The failure of the Pade resummation is not surprising because the Pade method assumes 
the approach to scaling according to (|2j). However, in the case of the Blasius equation 
the approach to scaling is as we will see in equation (p2|) in the next section. 



VII Variational Perturbation Theory for the Blasius 
Equation 

Variational perturbation theory for the Blasius equation fails to converge to the correct 
answer in the same way as for the instanton problem. We determined the leading strong- 
coupling term ( ^2|) up to 200th order and again it was impossible to find extrema, inflection 
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points, or higher derivatives that yield the correct resuh. Tables |lj and show the last 20 
strong-coupling coefficients 6g^'' and six orders of Richardson extrapolation. By determining 
the values of p and q we show why variational perturbation is likely to fail for this problem. 

Consider again the Taylor expansions for /„±i in (^4|) together with the Taylor series 
for fn-2 = fixn ~ 2a), namely 

fn-2 = f(xn) - 2/'(x„)a + 2/"(x„)a2 _ p"'{xn)a^ + lf""{xn)a^ ± ... . (76) 

Inserting these expressions into the difference equation for the Blasius problem ( [72| ) and 
translating back to the continuous function /(a;„) = /„, we get 



2e (^r'(x)a - lf""{x)a^ + + /(x) (^/"(x)a2 + \f""(x)a^ + = 0. 



(77) 



Next we transform back to the function y{x) — af{x) and assume the Taylor series 

y{x) = yo{x) + ayi{x) + a^y2{x) + ... . (78) 
To zeroth order in a we obtain 

2ey'^'{x)+yo{x)y'^{x)^0, (79) 



which is just the Blasius equation ([7 Of). The small parameter a, which is the lattice spacing, 
relates e and 5 hy a — \/e/5. Thus, if we evaluate the Taylor series (|7^) for the second 
derivative at the origin, we see that 

2/"(0) = y'm + ay'm + ... = ^-^^ + ^^y'm + ... . (80) 
Comparing this series to the original weak-coupling series 

y"(0) - (1 - 25 2^2 + ...), (81) 
we can now determine the leading power p/ q and the approach to scaling 2/q: 

= (5"^/2 (0.332O6 -I- (5"^/^eyi (0) -I- , (82) 



1 - 25 + 25^ 



so we obtain p = —2 and g = 4. 

Again we find that the approach to scaling 2/g — 1/2 lies just on the boundary of the 
open interval (1/2, 1), for which the proof of convergence ||l^ holds. This situation here is 
the opposite of the instanton case in that it sits at the lower boundary of the open interval 
in which variational perturbation theory works. 



VIII Large-order Behavior for the Blasius Equation 

The Blasius equation exhibits a large-order behavior which is a more subtle than for the 
instanton problem (p8|). The Blasius weak-coupling coefficients are not of Borel type; that 
is, the sign pattern is not alternating. Rather, the sign structure is governed by a cosine 
function with a frequency that is significantly different from tt. Remarkably, it turns out 
that a pure cosine cos(an) cannot reproduce all signs correctly. Up to 300th order the sign 
structure given by cos(an) is broken twice: The signs at n = 62 and at n = 212 are not 
correct if we optimize with respect to a. So we must consider an additional phase shift 
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cos(an + b). The parameter b turns out to be slightly smaller than tt, but it reproduces all 
300 signs correctly. 

In order to determine the numerical values of a and b we define 



fia,b) 



cos(an + b) ai^> 



N 

y . 

I cos(an + 6)1 |ai,„| 



(83) 



The sum ends at = 300 because this is as high as we can calculate using Maple; we know 
the first 300 weak-coupling coefficients ai j. For the correct values of a and b the function 
/(a, b) must be equal to 300. We then plot the function /(a, b) over the a-b plane and search 
for peaks. A careful study of the peaks yields values for a and b which allow the function 
/(a, b) to assume its maximum at 300. These numbers are given in Table |lj. 

The large-order behavior of the Blasius weak-coupling coefficients (unlike the large-order 
behavior of the instanton coefficients) has an additional overall factor cos(an -t- b), and we 
can now see that the remaining structure differs from the structure of the instanton weak- 
coupling coefficients. Dividing by the cosine, we observe that the coefficients 



•' cos{aj + b) 
irow factorially fast. Thus, we also divide hy j\: 



cos(aj -I- b)j\ 



(84) 



(85) 



The coefficients bj are unstable under a ratio test. That is, the ratio bj+i/bj decreases and 
then begins to oscillate. This is the inaccuracy that results from the delicate sign pattern 
of the first 300 coefhcients aij. 
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Tables 



j 




j 




1 


1 

8 


11 


2887747 
262144 


2 


1 

8 


12 


99392471 
4194304 


3 





13 


215798295 
4194304 


4 


11 

128 


14 


3781670831 
33554432 


5 


23 
128 


15 


8349041385 
33554432 


6 


295 
1024 


16 


1188129285795 
2147483648 


7 


589 
1024 


17 


2659104132291 
2147483648 


8 


39203 
32768 


18 


47890245452569 
17179869184 


9 


80723 
32786 


19 


108383753179167 
17179869184 


10 


1354949 
262144 


20 


39433620359113981 
274877906944 



Table 1: The first 20 weak-coupling coefficients aij for the instanton problem ( [l5| ) and 




N 


Sn 


N 


Sn 


1 


1 


11 


0.709998411 


2 


0.840896415 


12 


0.708235422 


3 


0.781934407 


13 


0.706789935 


4 


0.757237797 


14 


0.705659505 


5 


0.740759114 


15 


0.704734605 


6 


0.731210449 


16 


0.704006945 


7 


0.723927185 


17 


0.703419862 


8 


0.719045188 


18 


0.702964717 


9 


0.715146335 


19 


0.702610220 


10 


0.712308458 


20 


0.702349024 



Table 2: The first 20 Pade approximants for the solution to the instanton problem (19). 
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N 




N 


% 


180 


0.707530492 


190 


0.707471024 


181 


0.707524250 


191 


0.707465419 


182 


0.707518076 


192 


0.707459872 


183 


0.707511970 


193 


0.707454384 


184 


0.707505930 


194 


0.707448952 


185 


0.707499955 


195 


0.707443575 


186 


0.707494044 


196 


0.707438253 


187 


0.707488197 


197 


0.707432986 


188 


0.707482412 


198 


0.707427771 


189 


0.707476687 


199 


0.707422609 



Table 3: The last 20 variational strong-coupling coefficients feg^' from Eq. (^2[). 



order 


value for b^'^^ 


convergence 


1 


0.70640049 


decreasing 


2 


0.70639983200 


increasing 


3 


0.706399832082 


increasing 


4 


0.7063998320858658 


increasing 


5 


0.706399832085884411 


increasing 


6 


0.70639983208588446498 


increasing 



Table 4: Six orders of Richardson extrapolations for the strong-coupling coefficient 6g (fco) 
up to iV = 200 for the instanton problem. The last value is only 0.099% away from the 
correct limiting value 1/V2 = 0.7071067812... . 



order 


value for A 


convergence 


1 


-1.4998 


increasing 


2 


-1.500017 


decreasing 


3 


-1.5000011 


decreasing 


4 


-1.49999874 


increasing 


5 


-1.5000004 


decreasing 


6 


-1.499999893 


increasing 



Table 5: Six orders of Richardson extrapolations for the exponent A of the large-order 
instanton weak-coupling coefficients, based on the first 200 weak-coupling coefficients. The 
value A = —3/2 is quite plausible. 
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order 


value for K 


convergence 


1 


2.46692 


decreasing 


2 


2.4668283 


increasing 


3 


2.46682911 


decreasing 


4 


2.466829065 


decreasing 


5 


2.4668290597 


increasing 


6 


2.4668290635 


decreasing 



Table 6: Six orders of Richardson extrapolations for the inverse radius of convergence K 
of the large-order instanton weak-coupling coefhcients, based on the first 200 weak-coupling 
coefficients under the assumption that A = —3/2. 



order 


value for Bi 


convergence 


1 


0.0170837 


increasing 


2 


0.0170864 


increasing 


3 


0.017087 


increasing 


4 


0.0170893 


increasing 


5 


0.0170908 


increasing 


6 


0.0170922 


increasing 



Table 7: Six orders of Richardson extrapolations for the overall factor Bi of the large-order 
instanton weak-coupling coefficients, based on the first 200 weak-coupling coefficients under 
the assumption that K — 2.4482906 and A — —3/2. The value of Bi strongly depends 
on the numerical values for A and K. Changing K in the sixth decimal place influences 
the third significant figure of Bi. Also, all the Richardson extrapolations are increasing so, 
strictly speaking, we only have a lower boundary for Bi . Thus, the accuracy of Bi may not 
be very good. 



order 


value for B2 


convergence 


1 


0.119069 


increasing 


2 


0.119083 


increasing 


3 


0.119093 


increasing 


4 


0.119054095 


increasing 


5 


0.119054125 


increasing 


6 


0.119054146 


increasing 



Table 8: Six orders of Richardson extrapolations for the overall factor B2 of the large-order 
instanton weak-coupling coefficients based on the first 200 weak-coupling coefhcients and 
the same assumptions as in the case of Bi (see Table The value of B2 depends strongly 
on A and K. 
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j 




j 


«ij 


1 


-2 


11 


oUoDoDozooo 

5457375 


2 


2 


12 


6325029622 
637875 


3 


8 
3 


13 


487693745019181 
13408770375 


4 


-6 


14 


4774319527974167 
37819608750 


5 


184 
15 


15 


430321251088745734 
2212447111875 


6 


136 
9 


16 


796235344548876790517 
603998061541875 


7 


11062 
105 


17 


2249988054506764174584049 
6776858250499837500 


8 


8162 
225 


18 


178060537619150189817796 
14237097164915625 


9 


10557416 
14175 


19 


13224896152219729667498038639 
1301909768346024337500 


10 


57628622 


20 


121756993154067534451733120837029 


99225 


1153217968487557347375000 



Table 9: The first 20 weak-coupling coefficients for the Blasius recursion relation (Q) and 
([75|). Observe that the coefficients ai.j are not of Borel type (they do not alternate in sign). 
A cosine function with a frequency different from tt governs the sign pattern (see Sec. VIII). 



N 


Sn 


N 


Sn 


1 


0.5 


11 


0.3574632121 


2 


0.4204482076 


12 


0.3563326651 


3 


0.3948201830 


13 


0.3553848048 


4 


0.3819443732 


14 


0.3545795944 


5 


0.3742062309 


15 


0.3538882842 


6 


0.3690504811 


16 


0.3532891509 


7 


0.3653779673 


17 


0.3527655813 


8 


0.3626359060 


18 


0.3523046588 


9 


0.3605155915 


19 


0.3518961929 


10 


0.3588309707 


20 


0.3515320399 



Table 10: The first 20 Fade approximants for the solution to the Blasius equation ([70|). 
The sequence formed by the Sn converges extremely slowly. 



order 


value of y"(0) 


convergence 


1 


0.3445 


decreasing 


2 


0.3436 


decreasing 


3 


0.3430 


oscillating 



Table 11: Three orders of Richardson extrapolations for the Blasius equation (|7C|), based 
on the first 70 Fade approximants Sn- 
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N 


7 Af 

"o 


TV 


7 (N) 
% 


180 


0.33696017793094 


190 


0.33695971119646 


181 


0.33696012777085 


191 


0.33695966849139 


182 


0.33696007843082 


192 


0.33695962644843 


183 


0.33696002989308 


193 


0.33695958505396 


184 


0.33695998214034 


194 


0.33695954429471 


185 


0.33695993515575 


195 


0.33695950415774 


186 


0.33695988892292 


196 


0.33695946463046 


187 


0.33695984342591 


197 


0.33695942570058 


188 


0.33695979864918 


198 


0.33695938735612 


189 


0.33695975457760 


199 


0.33695934958540 



Table 12: The last 20 variational strong-coupling coefficients 6q for the Blasius equation. 
The very last coefficient is b^^^^ = 0.33695931237713, as opposed to the correct value 
y"(0) = 0.33206. 



order 


value for 6q^^ 


convergence 


1 


0.3369518 


increasing 


2 


0.336955563 


increasing 


3 


0.336955600539 


increasing 


4 


0.3369556008803 


increasing 


5 


0.336955600883462 


increasing 


6 


0.33695560088349232 


increasing 



Table 13: Six orders of Richardson extrapolations for the strong-coupling coefficient 
b(f^\ko) up to TV = 200 for the Blasius equation. The last value is 1.5% away from the 
correct limiting value t/"(0) = 0.33206. 



a 


b 


1.3941 


3.09 


1.3939 


3.11 


7.67830 


3.031 


7.67686 


3.130 



Table 14: Examples of the parameters a and b that give the first 300 signs of the Blasius 
weak-coupling coefficients correctly, assuming that the sign structure of the underlying large- 
order behavior is of the form cos(an -|- b). The last two values for a can be obtained 
approximately by summing 2n to the first two values. 
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Figure 1: The real part of the Pade approximants ^jv up to 200th order. Note that the 
approximants do not converge to the exact solution, which is represented by the horizontal 
solid line. The phases where the approximants become complex are marked by spikes. 
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Figure 2: The function from @ for N = 200 (solid line) and its second 

derivative with respect to ko (dotted line). The upper horizontal line equals l/-\/2, the 
correct limiting value of the instanton problem. All extrema of bg^^ are far from this value. 
Only the inflection point on the right-hand side comes close. The value for fco, for which the 
second derivative vanishes, is fco = 18.42510. Substituting that number into the function 
&g^-'(fcp^-'), we obtain in the 200th order b'^^^^'' — 0.707417. The corresponding Richardson 
extrapolations can be found in Tabic 0. 



